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The probability that a point is to one side of a curve in Schramm- 
Loewner evolution (SLE re ) can be obtained alternatively using boundary 
conformal field theory (BCFT). We extend the BCFT approach to treat 
^ , two curves, forming, for example, the left and right boundaries of a cluster. 

■ This proves to correspond to a generalisation to SLE(«;,p), with p = 2. 
' We derive the probabilities that a given point lies between two curves or to 

f^*) , one side of both. We find analytic solutions for the cases k = 0, 2, 4, 8/3, 8. 

0^ ' The result for k = 6 leads to predictions for the current distribution at 

| the plateau transition in the semiclassical approximation to the quantum 

*0 ! Hall effect. 

O 

^ 1 Introduction 

A large class of two-dimensional lattice models may be described in terms of 
^ a gas of non-intersecting loops. Examples include the boundaries of Ising spin 

JL" ■ clusters and percolation clusters, the boundaries of the clusters in the random 

. i-H ! cluster representation of the Potts model, the level lines of solid-on-solid models 

■ of the roughening transition, as well as dilute self-avoiding walks and polygons. 
These are in fact all special cases of the lattice O(n) model. 

The continuum limit of these curves is conjectured to be both scale and confor- 
mally invariant (a statement which may be made precise and which has been 
proved in a few cases.) A description of the whole ensemble of such curves is dif- 
ficult, and instead it is simpler to focus initially on a single curve connecting two 
given points on the boundary of a simple connected domain, whose existence is 
guaranteed by the boundary conditions. The conformally invariant measure on 
such curves is conjectured to be given by Schramm(stochastic)-Loewner Evo- 
lution (SLE) PP with parameter k, where n = — 2cos(47t/k) with 2 < k < 8. 
Many of the previously conjectured scaling dimensions of the O(n) model, as 
well as other properties such as crossing formulae, have been derived using SLE 
0E1 El- The stochastic description naturally leads to second-order linear differ- 
ential operators. Some of the scaling dimensions and correlation functions are 
given in terms of suitable eigenvalues and eigenfunctions of these. 
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An alternative description of the continuum limit of the critical 0(n) model is in 
terms of conformal field theory (CFT). This focusses on correlation functions of 
local scaling operators. Within CFT there is a correspondence between such op- 
erators and states in the radially quantised theory. These may be organised into 
irreducible representations of the Virasoro algebra satisfied by the generators 
L n of infinitesimal conformal transformations. 

A connection between these two pictures was made in Refs. [SJ 0. As was 
conjectured in 19840, conditioning the CFT partition function on the existence 
of a curve starting at a given boundary point is equivalent to the insertion of 
a boundary operator 1 which corresponds to a Virasoro representation with 
a level two null state: i— 2 1<^2) = otL^i^fa) ■ As was also shown in 1984|S], this 
implies that the correlators of this operator satisfy certain second-order linear 
differential equations. These are the same as those coming from SLE K , with the 
identification a — k/4. 

However, CFT describes not only single curves but also many. Indeed the 
correlator 



conditions the partition function on the existence of N such non-intersecting 
curves, hitting the boundary at the points {xj}. Such a correlation function 
satisfies 2N linear BPZ equations, one for each xj. Note that the order in 
which they link up is not specified in the above. It may be made more precise 
by assigning O(n) labels to the t^s: only operators with the same label can then 
be connected by a curve. Alternatively, as we shall discuss further, the different 
ways the curves link up correspond to imposing different boundary conditions 
on the differential equations. 

It is therefore straightforward in principle to derive within CFT many different 
results relating to how N such curves cross a given domain. In practice, this 
becomes technically prohibitive. In this paper we consider the simplest non- 
trivial case with N = 2. An interesting application of this for k — 6 (percolation) 
is the following problem, illustrated in FigQ consider critical site percolation 
in the upper half plane, and suppose that all the sites on the real axis are 
constrained to be white, except that at the origin, which is black. Moreover 
this site is conditioned to be connected to infinity by black sites, that is, it is 
part of an incipient infinite black cluster. The boundaries of this cluster then 
define two curves of the type we are considering. 

One of the simplest SLE results for a single curve, due to Schramm^]], gives 
the probability that a given point £ in the domain lies to the left (or right) of a 
single curve. While this had not in fact previously been discussed within CFT, 
its derivation from this point of view is st r ai ght forward [HI] : the probability is 
given by the ratio of conditional expectation values 



where $(£) is the indicator function that the curve passes to the left (right) of 
£. [This behaves to all intents and purposes like a local operator with scaling 

x We use the notation <j> n rather than the standard <j> n ,i or 4>i, n because the paramctrisation 
in terms of re does not distinguish between the last two. 
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Figure 1: An incipient infinite cluster containing a point on the boundary in 
critical site percolation on the triangular lattice. The scaling limit of the two 
curves forming its boundary is described by the results of this paper. 



dimension zero.] Specialising (without loss of generality, because of conformal 
invariance) to the case when the domain is the upper half plane and x 2 — > oo, 
the denominator becomes trivial and the numerator satisfies a BPZ equation 
with respect to x\. 

The generalisation of Eq. J2J to N curves is straightforward. For example, for 
two curves in the upper half plane conditioned to connect (x\, oo) and (x 2 , oo) 

p(f] = ($(£)</>2 (xi )jh {x 2 )(f>2 M02 (oo)) ,„s 

where now $(£) is any of the indicator functions that £ lies to the left, or right, 
of both curves, or between them. Both numerator and denominator satisfy 2nd- 
order BPZ equations with respect to x\ and x 2 : the boundary conditions pick 
out which case is being computed. 

These coupled partial differential equations are already too difficult to solve in 
closed form. However, they simplify in the limit when S = x 2 — x\ <C |£|, by 
a well-known property of CFT called the fusion rules. These state, roughly 
speaking, that in this limit the operator product ^2(^1)^2(^2) may be written 

<M*i) ■ 02(^2) = s^m^i) + s a2 Mxi) , (4) 

where the correlators of 4>\ and 03 satisfy respectively first and third order 
equations. In fact c/>3 corresponds to a Virasoro representation with a null state 
at level 3. The null state condition may be computed explicitly, and hence 
the third-order equation. (This could also be obtained directly from the two 
2nd order equations.) We also argue that if the curves are conditioned to go 
to infinity (rather than there being a single curve connecting x\ and x 2 ) this 
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picks out the second term in Eq. fllf. The resulting correlation functions now 
depend only on arg(£— x\) and may be written as indefinite integrals of ordinary 
hypergeometric functions. 

The layout of this paper is as follows. In the next section we recall the derivation 
of Schramm's formula for a single curve using both SLE and CFT. SecGDcontains 
the main body of the CFT calculation for two curves. In Sec.0]we present the 
results, in graphical form for n = 6, 8/3. For some other values of k they may 
be expressed in terms of elementary functions. The limiting cases k = and 
K = 8 are interesting. 

In Sec. we apply the results for n = 6 to the semi-classical limit of the quan- 
tum Hall plateau transition, where electrons in a strong magnetic field move in 
random scalar potential. If Coulomb interactions may be neglected, the guiding 
centres of the electrons approximately follow the level lines of the potential, 
which, in the scaling limit, are percolation cluster boundaries. At the criti- 
cal point the half-plane geometry may be conformally transformed into a long 
strip, and our calculation gives information about the paths followed by the 
conduction electrons, and hence the mean current density across the strip. 

Finally in Sec. we discuss the equivalence between the BPZ equations of CFT 
and a postulated generalisation of SLE to multiple curves. The growth of a single 
curve, conditioned on the existence of the others, is described by a variant of 
SLE known as SLE(k, p) (with p = 2), as asserted in several recent studies of 
multiple curves [TT1U21E3- The CFT equations also suggest [HE| that [t is 
possible to describe the joint measure on all the curves by a 'multiple SLE' in 
which they are grown simultaneously, as shown precisely in |12j . 



2 Schramm's formula 

In this section we review the computation of the probability that a given point 
£ lies to the left(right) of a single curve, from the points of view of both SLE 
and CFT. 



2.1 CFT method 

As discussed in Sec. ^ the probability is given as the ratio of correlators in 
Eq. 10). We take Xi — > x and X2 large (and eventually to infinity). Consider the 
effect of the infinitesimal conformal transformation z z + e/(z — x), which is 
implemented by inserting into each correlator a factor 

eT(z) dz (s) 



c 



z 



x 2ni ' 



where C is a contour surrounding x (but not £ or x%), together with its re- 
flection in the real axis. This may be evaluated in two ways: by shrinking the 
contour around z = x and using the fact that the 0((z — x)°) term in the 
operator product expansion of T{z) with 4>2{x) is (by definition) L_i<i>i{x) = 
(k/ A)L 2 _ 1 4>2(x) = (k/4)9^</>2(x); or by wrapping the contour around £ and Xi- 
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The effect on $(£) is just to shift £, while the effect on (^2(^2) is negligible as 
x-2 — * 00. 

Equating these two ways of evaluating the insertion gives, as X2 — > 00, the BPZ 
equation 

'k 9 2 _ r 1 d 



4&» l2RC 



P(£,x)=0, (6) 



where P(£, x) is either probability. Since these only depend on the angle which 
£ — a; makes with the axis, or cquivalcntly the variable t = Re(£ — x)/Im£, this 
partial differential equation reduces to an ordinary one, of Riemann type. The 
boundary conditions are Pi e ft — > as t — > +00 and Pi e ft — ► 1 as i — » —00. The 
equation itself allows these two possible asymptotics - an explicit calculation 
determines the exponents as \t\ — > 00 to be |i|~ 7 with 7 = or 7 = X2 = 
(8 — At)/ k. The latter is the boundary 2-leg exponent: it arises as t — > +cx) since 
the point £ traps the curve against the real axis so that, in effect, two mutually 
avoiding curves of the O(n) model emerge from that point. 

These boundary conditions then fix the solution to be 
p _ 1 r (!) , F A 4 3 

Plcft - 2 ~ 0?r(^) 2 Fl( 2'«'2'~* } ' (7) 

Note that, because Pi e ft + -Fright = 1 is a solution to the equation, any other, 
including Eq. J7J , may be written as a quadrature of an elementary function. 



2.2 SLE method 

We summarise the theoretical physicist's version|2 E] of Schramm's original 
argument. In SLE, the curve in the upper half plane from a point ao on the real 
axis to infinity is considered as being grown dynamically, introducing a fictitious 
time variable t (obviously distinct from the variable t defined above.) Let K t 
be the set consisting of the curve as grown up to time t (as well as all points 
enclosed by the curve and between the curve and the real axis) so that the 
complement of this set in the upper half plane is simply connected. Let gt{z) 
be the (unique) conformal mapping of this complement to the whole upper half 
plane, normalised such that gt(z) = z + 0(1/ z) as z — > 00. The coefficient of 
the 0(1/ z) term is increasing with t, so 'time' can be reparametrised so this 
coefficient is exactly 2t. The image of the growing tip of the curve under g t is a 
point at on the real axis. Loewner showed that the time-evolution of gt satisfies 

d 9t(z) = 2 
dt gt(z)-a t ' 1 ' 

Any suitably continuous function at generates a curve. SchrammpQ showed that 
if this process is to generate a conformally invariant measure on curves, the only 
possibility is at = -\J~RBt + ao, with B t being a standard Brownian motion. 

Now consider the problem at hand, with a curve 7 connecting ao to 00, and 
a given point £ away from ao. Evolve the SLE for an infinitesimal time dt. 
The function g^t will erase a short initial segment, and map the remainder of 
7 into its image 7', which, however, by conformal invariancc, will have the 
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same measure as SLE started from a^t = ao + \ /f KdB t . At the same time, 
£ -> £' = £ + 2eft/(£ - a ). Moreover, 7' lies to the left(right) of £' iff 7 lies to 
the left(right) of £. Therefore 

P(£; ao) = (P(£ + 2tfc/(£ - ao), «o + v^t)} , (9) 

where the average (. . .) is over all realisations of Brownian motion dB t up to 
time dt. Taylor expanding, using (dBt) = and ((dB t ) 2 ) = dt, and equating 
the coefficient of dt to zero then gives exactly the CFT Eq. © if we set ao — > £. 



3 The differential equation for two curves 

According to CFT the probability that a given point £ lies to the left, between 
or to the right of two curves starting at points x\ , X2 on the real axis is given by 
a ratio of correlators as in Eq. J^J. The numerator and denominator each satisfy 
BPZ equations with respect to both x\ and X2- For general values, these are 
discussed further in Sec. |S| However, explicit analytic progress is only feasible 
in the limit when 6 = x-i — x\ — > 0, and we now treat this using established 
properties of CFT. 

The operator product expansion of two <f>2 operators is constrained by the fusion 
rules to have the form in Eq. (0J). This means that every solution of the coupled 
second order BPZ equations may be written in this limit as a linear combination 
of functions 6 aj Fj(£, x, 5), with j = 1, 3 and x = (xi + X2)/2, and the functions 
Fj having a regular power series expansion in 5. The values of the ctj are 
determined by the differential equations to be a\ = — 2/i2 = —(6 — k)/k and 
«3 = hs — 2h 2 = 2/k. 

In general, the dominant behaviour as i5 — > is given by a±, corresponding to 
the first term </>i in the OPE, Eq. J3J. This is just the identity operator, and it 
is straightforward to see that the corresponding solution F\ is in fact a constant. 
The physical interpretation of this is that Xx,x% are overwhelming likely to be 
the end-points of the same curve as 5 — > 0, which makes a very small excursion 
into the upper half plane. It has no effect on conditional probabilities of events 
further away. In order to condition the two curves each to go to infinity, we 
must therefore impose the condition that this term is absent in the solution. 

This leaves the term in the OPE coupling to the ^3 operator. The leading 
behaviour of the probability function, Eq. @, in this limit is therefore given by 
the ratio 

P = km mpmm . (10) 

In the limit y — > 00 the denominator is trivial, going as y~ 2h3 and serving only 
to make the whole expression finite. According to CFT [HlEjj the correlator of 
the (f>3 operator in the numerator satisfies a third-order equation of the form 



2Rc 



(£ - x) 2 d£ 



2^Rc 



1 



d <9 3 \ 



where fi and A arc defined by the level 3 null state condition, 
(L_ 3 + /ui-ai-i + AL 3 ^)!^) = 0. 
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In Appendix 1X1 it is shown that 

^ = -T' A = ; 7TT) v ( 12 ) 
h 3 h 3 (l + h 3 ) 

where h 3 is the conformal scaling dimension of (j> 3 . A derivation of Eq. 1)11(1 is 
included in Appendix^] As with the one curve case, the function P is expected 
to depend only on its angle from the imaginary axis, or equivalently, the variable 
t = (u — x)/v, with u and v coming from £ = u + iv. This can be used to write 
the partial differential equation as the following ordinary differential equation 
in t 

d 3 P 2 M t d 2 P (3- M )t 2 -(l + ^) dP 

dt 3 t 2 + 1 df 2 (t 2 + 1) 2 dt 1 ' 

After some algebra this may be further rewritten in terms of Q(t) = dP(t)/dt 
and the variable s = —t 2 as 

^d 2 Q A-(A-2/i) 5 dg (3-/i).s + (l + M ) ^ 

s(l - s)-— H — 1 — — Q = , 

ds 2 2A ds 4A(1 - s) 

which is of Riemann form, see Chapter 4 of |17j . with the following exponents: 

• as s « 0, Q - s°, s 1 ' 2 , 

• as s ps 1, Q ~ (s - l) 1 " 8 ^, (s - l)- 8 / K , 



• as s 



CO, Q ~ S- 4 / K , s -(24-k)/2 K 



The solutions to the ordinary differential equation therefore take the form 

( 1+i 2)*-i ■ (W) 



3.1 Boundary conditions and solutions 

Recall from Section that there are three possibilities for the position of a point 
in the upper half plane: it may be to the left of both curves, between them or 
to the right of them. Each case corresponds to different boundary conditions on 
P. Consider firstly the case that the point lies to the left of both curves. The 
boundary conditions for this case are 

lim flrft(f) = l, 

t— > — OO 

lim P lcit (t) = 0(t-W K -V). 

t— >+oo 

The first condition comes from insisting that the probability of a point on the 
negative real axis being to the left of both curves is one. The second condition 
is best understood with reference to Figure below, from which it can be seen 
that the limit t — -> oo should have the exponent corresponding to a four-leg 
operator on the real axis, namely 24 = /15 = 24/k — 2. 
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Figure 2: Pi e f t as t — > oo 



Recall from the argument below Eq. H13[). that 

^ c 

with Q(t) given in Eq. (|14|) and c a constant, to be determined by the boundary 
conditions. In order to apply the boundary conditions at large t, the hyperge- 
ometric functions must be analytically continued for \t\ > 1, see for example 
15.3.8 in ^Hj- The expression for P then takes the form 



Pit) = A 



iXi)r(l-f) 



r(i-i)r(-|)(i + u2) ^ 
r(|)r(f-|) i 



T 2 



4 4 

k' k 



r(^l)r(M)(i + « 2 )i 
r(|)r(i-f) 



r(§-£)r(i-£)(i + ti»)* 
r(|)r(f-i) 



2^1 (1 - 
2^1 (1 



4^ 

K 

4 4 1 

K K 2 



1 8 
2' « " 

; l 



2' 1 + m 2 

? 1 - 



2' 1 



r) 



r(i + i)r(i) ( i + u 2)f 



Fi(-. 



4 1 

k'2 



k' 1 + u 2 ' 

4 3 8 1 
k'2 _ k'TT" 



r) 



dit . 



The first and third terms in this large t expansion approach zero as £ 2 - 24 / K ; while 
the second and fourth term approach zero as £ 1 ~ 8 / K . For the range < n < 8, 
which is the range of physical interest, the first and third terms fall off more 
quickly with increasing t. In order to satisfy the second boundary condition, the 
coefficients of the dominant second and fourth terms must ensure cancellation 
as t — > oo. This uniquely determines the constant B as 



B = -2- 



r(i + t)i\i) 
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The solution follows immediately, 



PeftW 



f t °°S(t')dt? 



where 



S(t) 



k' 2' J r(|+^)r(-|+f) 



2£(i±lHi) / p n 



4 3 _ 4. 3. 

re> 2 re ' 2 ' 



+ 

It may be simplified using the identity 

//N , 2 2 - 8 / K 7rr(^ - 1) 
S(t')dt' = , N " . 

r(±)r(£) 

see for example Eq. (7.512.10) in JHj- The solution then takes the more elegant 
form 

r(!)r(f) 



Plcft(t) 



2 2-8/ K7rr (il _ !J 



S(i')di' . 



It is a simple matter to derive P r i g ht since P r i g ht(i) = Pieft( — thus 

r(!)r(t) 



Pright W = 



2 2-8/ K7rr (il _ !) 



5(t')dt' . 



(15) 



(16) 



The remaining solution is P m iddic- This can be derived in two ways. Firstly, for 
the total probability to be unity, Pi e f t + P m iddic + fright = 1- Subtracting the 
two previous solutions from one leads to 



P 



middle 



(*) = i - 



2 2-8/ K7rr (i^ _ !) 



S(t')dt' 



S(t')dt' 



(17) 



Alternatively, the solution can be written as the integral of the odd term in 
S(t) only, since we expect the solution to P m iddie to be an even function of t. 
Applying the boundary condition P m iddie(i — > oo) = determines P m iddie up to 
a multiplicative constant 



Pmiddlo(£) — D 



4 3 _ 4. 3. 

2 k ' 2 ' 



-u 2 )du 



(18) 



The constant D may be found by equating this expression with Equation l|17l) 
at t = 

r(£)r(£) 



D = 1 — 



2 i-8/« 7r r(if - 1) 7 



5(i')di' . 
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Below, in Figure © are the results for k = 8/3, conjecturally corresponding to 
the scaling limit of two mutually avoiding self-avoiding walks starting near the 
origin, and k = 6, corresponding to the scaling limit of the percolation problem. 
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4 Special Cases 

In the special cases k = 0, 2, 8/3, 4, 8, explicit solutions may be derived for Pi e ft, 

P middle and P r ight- 



4.1 k = 

With k = 0, the curves are deterministic straight lines. They start at the origin 
and proceed at an angle of 7r/3 radians from each other and from the real axis, 
as in Figure 10}. Actually this is a rather singular limit of the equations, and the 
above result is more easily understood following the multiple SLE interpretation 
of Sec. El 




Figure 4: k = 



4.2 K = 2 

The solutions are written in terms of the function S(t). In the case of k = 2, 
this takes the form 



s(t) = u 



1 



(1+t 2 ) 3 
1 

(1 + * 2 ) 3 



l + 5t 2 



32< /i(13 + 15t 2 ) + 3(1 + 6t 2 + 5t 4 ) arctan(t)~ 



3tt V 



16(i + t 3 ) 
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Inserting this expression into Eq. I|15|) . (|16fl and (|17|l and performing the inte- 
grations, the following solutions are obtained 

1 1 



4 97T 2 (l+f 2 ) 3 



(-16 - 9t 2 + 9t 4 ) - 9ir(t 3 + t 5 ) 



+ 9(1 + t 2 ) arctan(t) [2t 3 - tt(1 + t 2 ) 2 + (1 + t 2 ) 2 arctan(t) 



(19) 



1 



bright (*) = T 



1 



4 9vT 2 (l+f 2 ) 3 



(-16 - 9t 2 + 9i 4 ) + 9tt(F + i") 



+ 9(1 + r) arctan(t) 2i 3 + tt(1 + t 2 ) 2 + (1 + r) 2 arctan(i) 



-Pmiddlo(i) — 7 — 



(-16-9i 2 + 9i 4 ) 



1 2 

2 ~ 9tt 2 (1 +< 2 ) 3 

9(1 + t 2 ) arctan(t) (2i 3 + (1 + i 2 ) 2 arctan(t) 



(20) 



(21) 



4.3 « = 8/3 

Analytic solutions may be obtained for the special case of the scaling limit of 
two self-avoiding walks, conjectured to correspond to k = 8/3. In this case, the 
function S(t) is given by 



1 



(1 + i 2 ) 2 



1 1 9N 3nt 



S(t) = n , , 2 , 2 2 Fi(2,--;-;-^)- 



2'2' ' 4 



1 / 1 3t ^ 3 3?rf\ 

— arctan(t) + - - — 



(1 + t 2 ) 2 V 2(1 + 1 2 ) 2 w 2 4 / 

The solution for Pi c f t follows as 

2 2-8/ K7r r(f - 1) J t 

16 Z* 00 1 / 1 3*' , 3 3nt' 

vto] v (iT^(-27TT^) + T arctan(0 + 2- — 

16 r£'(13 + 15i /2 ) + 3(l + 6i /2 + 5£' 4 )arctan(f) 3tt 1 



*W*) = ^_ s/ :_w!2 ^ / S(i')df 



15tt 



16(1 + i' 2 ) 2 8 1 + i' 2 



1 16 i(13 + 15i 2 ) + 3(1 + 6i 2 + 5t 4 ) arctan(i) 2 1 



2 15tt 16(1 + t 2 ) 2 51 + t 2 ' 

This may be simplified to 

— 2£(13 + 15£ 2 ) + 3tt(1 + 6t 2 + 5t 4 ) - 6(1 + 6t 2 + 5t 4 ) arctan(t) 

loftW ~ 30tt(1 + t 2 ) 2 ■ 

(22) 

The other solutions are 

2t(13 + 15i 2 ) + 3?r(l + 6i 2 + 5t 4 ) + 6(1 + 6i 2 + 5t 4 ) arctan(i) 

p «Mt) somTTW ( } 

p *^® = 50Tey (24) 
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4.4 k = 4 



With k = 4, takes the form 



1 



1 + t 2 



3 „ 1 2, n „ ,„ 1 3 o, 

■■^ I D+ _ IT 1 . /'O +2^ 



= 2F 1 {-,0,-,-t 2 )+Bt 2 F 1 [2,-,-,-t 



2' 1 2' ' v 2' 2' 



where £? = —Ajir. Substituting for the hypcrgcomctric functions in terms of 
elementary functions leads to 



S(t) = 



1 + t 2 



^ 2t ( 1 arctan(t) 



7T \l + t 2 



Inserting this form into the expression for Pi e f t 

1 1 arctan(t) arctan 2 (i) , . 

= 1 - ^iT^ — ir 1 + ■ (25) 



The other solutions arc 



1 1 arctan(i) arctan 2 (i) 



1 2 2arctan 2 (i) 

-Pmiddlc(i) - 7T + 2C1 i +2\ 2 ' 

2 TT Z (1 + t £ ) 7T 2 

In terms of the angle <p from the imaginary axis to the point £, which is to say 
(f> = arctan(u/w) = arctan(i), 



-Pright 0) 
-P m iddle(<A) 



l cos 2 (0) 4> 4> 2 

4 7T 2 7T 7T 2 

1 cos 2 (0) (j) 4> 2 

4 7T 2 7T 7T 2 

1 2cos 2 (» 20 2 

2 + ^ ^ 2 "' 



4.5 k = 8 

The case k = 8 is subtle and requires careful treatment. For the single curve, 
taking the limit as K — » 8— at fixed t yields an expression which fails to satisfy 
the boundary conditions (Hj- Instead, the probability of any point not on the 
real axis being to one side of the curve is everywhere equal to a half. This 
however makes physical sense since for k = 8 the curve is space-filling. 

For the case of two curves there is a similar boundary condition violating solution 
which is obtained by continuing the general solution for finite t to K = 8. There is 
a second solution, however, which follows from solving the differential equations 
at k = 8 and which does satisfy the boundary conditions. 



First, consider the general solution in the limit k — > 8. Using the usual definition 

2r(i+-|)r(A) 



jp (1 , 4 1 4.1. .2\ 2r ( 1 +-) r (~) 4 TP n i 4 3 4.3. ±2\ 
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the solution is 



fleft(*) = 



22-8/ K7rr (^ _ !) 



In anticipation of taking the limit n —> 8, define k = 
small. Rewriting the above expressions in terms of e 



S(t')dt' . 

= 8 — e and assume that e is 



S(t) = 



fleft(*) 



2^1(1 



16' 2 



_e_. 1. 

16' 2' 



g£(l±AHi±A) / p r3 
r ( 1 +Ts) r (iij) t2 ^ l( -2 
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,1 



_e_. 3. 

16 ' 2 ' 



r (! + T6-)r(i + f) 



(1 + 1 2 )* 
5(t / )dt / . 



2i-* 7r r(i + f§) 

The two hypergeometric functions in S'(i) can be analytically continued for i > 



5(t) 



1 



" r(| - £)r(- 



e e 3 c 
2Fl(1+ 16'16 ; 2 + ° 



1 



r(i + ^)r(^) (1 + 



e 

16' 



e 1 
16' 2 



8' 1 4 

e 1 

8' 1 + i' 2 



t' 2 



_e_ 1 , _e_. 3 

16' 2 "r" 16' 2 



l+t' : 



r(i + ^TO)r(i-^)r(-^) 
. r(| + |) t' ^(1-^,-^1 



(l+f' 2 ) 1+ TF 



Each of these four terms are to be integrated from t' = t to 00. This can be 
done by expanding the hypergeometric functions and integrating term by term. 
Consider the integral of the first function in S(t) for t > 



i) 



r(!-^)rH 

V5rT(-i 



1 



f) 



16 > 



(l+t' 2 ) 1+ T§ 
1 



2^1(1 



OO 



e e 3 
16' 16' 2 + 

(1 + T6-XT6-) 



1 



1 + t' 2 



)dt' 



16 ' 



_ ,1 3e 3 ,. 

= -t 2 F 1 (-,l-\ ;-;-i 2 ) 

y 2 16 2 ; 



(l + i' 2 ) 1+ ft 

OO 

(l + i' 2 )^ 

0?r(i + ff) 



At' 



2r(i + i) 



(| + |)1! 1 + t 2 
ef(t,e), 



dt' 



where /(e, i) is a convergent function of its arguments for all t in the limit e — > 0. 

After integration, the third term gives e 2 g(e, t), where g(e, t) is also a convergent 
function of its arguments for all t in the limit e — ► 0. 

The integral of the second and fourth terms should be considered together. They 
contribute 



v^n\ + 1) 

r(i + £)r(£ 



1 



(l + i' 2 )Hra (1 + <' 2 ) 1+ ts 



dt' + eh(e,t), 
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where h(e,t) is a convergent function for all t in the limit e — ► 0. The integral 
of 5(t) is therefore 

f S(t')df = -* 2j P 1 (|,l + f|;|; -t 2 ) + + e/ ^' £ ) + ^(t, e) 



v^a + 1) 



di' + e/i(e,t) 



r(i + ^)r(^) y t L(i + t«)i+A (1 + 

Taking the limit e — > 0, the first two terms become 

— arctan(t) + — 

and the terms involving f(e, t), g(e, t), h(e, t) have coefficients which tend to the 
zero. The remaining term is 



lim 



i a 



~o r(i + ^)r(^) J t L(i + t «)H* (i + i' 2 ) 1+ ^ 



where j(e, i) is a convergent function of its arguments for all t in the limit e — > 0. 
The integral is finite, so the expression goes to zero as e — > 0. 



Using 



lim 



r(| + ^)r(i 



e^o 2 1 -i 7 rr(i + f§) 2tt' 

allows the solution for Pi e a{t) and related expressions to be deduced in terms 
of the angle <f> defined in the previous subsection, 

fl B ftW = (28) 

bright (0) = J + ^ , (29) 

P mi ddic(0) = ^ • (30) 

Note that these are valid only for <j> ^ ±ir/2. The limits K — > 8 and |i| — > oo do 
not commute. 

Now take k = 8 from the beginning. The differential equation has solutions of 
the form 

fl eft (t) = ^arctan(t) + B . 

If A and -B are chosen to be —1/ir and 1/2 respectively, this satisfies the bound- 
ary conditions. Then in terms of the angle <f>, 

PleftW = ^ - - , (31) 
PrightW = £ + - , (32) 

P,„iddio(0) = . (33) 

This is not the analytic continuation of the first solution to k = 8. Mathemati- 
cally, this may be traced to the fact that the limits k — > 8— and t — > +oo do not 
commute, yet we have to impose the boundary condition on Pi e f t at t = +oo. 



15 



Physically, these two solutions appear to lead to different pictures. In the first 
case, the probability that a given point lies between the two curves is exactly ^. 
This may be understood 2 in terms of the physical picture of the curves being 
the boundaries of two disjoint uniform spanning clusters which are separated 
by some random simple curve 7 SO p- Each curve fills the entire region to the 
left (right) of 7 scp . Within each region, the probability that a given point lies to 
the left (right) of the curve is \ as before, but the probability that it lies in this 
region, to the left(right) of the scparatrix is |(1 =F 20/7r). Thus the probability 
that a point lies to the left of both curves is \ - ^(l~2<fi/ir), while the probability 
this it lies between them is 

| - | (1 - + I - I (1 + 20/tt) = i . (34) 

In the second solution, the probability that a point is between the curves is zero, 
hence the area between them vanishes. Physically, this could correspond to two 
mutually avoiding dense polymers, which also correspond to k = 8. In this case, 
however, they find it entropically favourable closely to follow each other. The 
probability that a point lies to one side of this composite curve is the same as 
for a single curve with k = 4, see ||]. This may be understood in terms of 
the stochastic interpretation of Sec. El For k = 8 the points x\ and X2 almost 
certainly collide in finite time, after which it is necessary to prescribe how to 
continue the solution. One possibility is that the two points coalesce, in which 
case their centre of mass describes a Brownian motion with k' = \n — 4, as in 
the second solution. Another possibility is that they are conditioned never to 
collide, which then presumably corresponds to the first case. 



5 Application to the quantum Hall transition 

The Quantum Hall effect is observed in two dimensional electron gases in semi- 
conductors with magnetic fields applied normal to the plane. Donor ions are 
spatially separated from the electron gas in order to increase the mobility of 
the electrons in the sample. The ions are positively charged and the resulting 
Coulomb potential in the electron gas can be modelled as a random potential, 
V. A semi-classical approximation may be used in the limits of slowly varying 
potential on the scale of the magnetic length I = ^/TicjeB and strong magnetic 
field, defined as 

TieB 

nuj c = >C e_F , 

mc 

where is the Fermi energy of the system. In this limit, the eigenfunctions of 
electrons are large only around constant energy surfaces of the potential, V, see 
|21j . The eigenfunctions can be approximated by 

i>{u 1 v) = C{u) Xn {v)e^\ 

where u is the length along the constant energy surface and v is the distance 
normal to it. C is a normalisation factor and Xn is the nth harmonic oscillator 

2 We are grateful to W. Werner for pointing this out. 
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function 

Xn = H n (-)e il 7 . 

Choose the zero of the random potential to be the spatial average of the potential 
and assume that (V(r)V(r + 5r)) goes to zero for 5r b. The requirement that 
the potential varies slowly compared to the magnetic length is equivalent to 
b >> I. If all points where the potential is greater than a value E are coloured 
white and points where V < E are coloured black, the lines of constant potential 
will be given by the boundary between the two coloured regions. This model 
is believed to be in the same universality class as lattice percolation. Thus, 
electrons will move along the boundaries of percolation clusters. In general, 
the boundaries will form closed loops, but for the critical value E = 0, their 
mean size diverges and they will be locally described by SLE with k = 6. The 
appropriate boundary conditions for the percolation picture are that the top 
and bottom edges should be coloured white since the potential is effectively 
infinite there. 

As a function of y, the distance across the sample, the energy of the Landau 
levels follow the form of the potential and those states with E < /i, the chemical 
potential, are occupied, as shown in Figure JSt). Diamagnetic currents flow, 
both around the closed loops and along the extended cluster boundaries with 
E = 0, as shown in Figure On connection of the leads to the ends and 

the application of a potential difference, the current distribution will change 
from that of the equilibrium case. The currents flowing along the extended 
cluster boundaries will be affected the most, since these are the only paths 
which can carry a net current along the length of the sample without potential 
tunnelling. The Fermi energy of electrons moving along the bottom cluster will 
be the equilibrium level plus eU/2, since these electrons flow from the electron 
reservoir at the negative battery contact. The Fermi energy of the electrons in 
the upper cluster boundary will be that of the equilibrium case minus eU/2, since 
these electrons flow from the electron reservoir at the positive battery contact. 
Relative to the equilibrium case, therefore, an extra current, 51 flows along the 
lower cluster boundary and, to a first approximation, 51 less current flows along 
the upper cluster boundary. The averaged change in the current distribution, 
compared to the equilibrium case, will be given by the spatial average of the 
extended cluster boundaries in the percolation picture (k = 6). We may relate 
the strip geometry of the Hall bar experiment to the half-plane discussed in this 
paper by the conformal mapping z — > z' = w(z) with 

w(z) — — ln(z) , 

which maps the upper half plane to an infinite strip of width L. Lines with 
constant angle, parametrised by their value of t = u/v, are mapped to lines of 
constant distance, y, from the bottom of the infinite strip, given by 

L / 1 \ 

y = — arctan(— j . 

7T t 

A little thought then shows that the mean extra current density flowing along the 
upper boundary curve at height y in the sample is proportional to the derivative 



17 




Figure 5: The spatial distribution of cigenstates, showing occupation before and 
after a potential difference is applied between the ends of the sample. 



with respect to y of the probability that y is above both curves. Similarly, the 
extra current density flowing along the lower curve is — dPbeiow/dy. Thus, 



dPleft dP ri 



ght 



X 



dy dy 

22-s/ K7r r(f - 1) l 



tan(f)^ 



tan(^) 



tan(7r-^y 



where, as previously, 

_ 1 . 1 . +2\ 



p f 1 I 4 i 4.1. .2\ 2r ( 1 + 7:) r (-) 4 C i 4 3 4.3. ±2\ 
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(b) 




Figure 6: The top diagram depicts the diamagnetic currents. The figure below 
shows the Hall bar with leads attached and a voltage applied. The arrows in 
both diagrams represent the direction of electron flow. 



For as = 6 this becomes 



IOC 



(1 



) S(- 



tanm^L Han(^) 



-) + S( 



tan(7r 



L ) 



with 



S(t) 



3' 2' 



2r(j)r(|) 



t 2-F\(§ , 



3. 
2 ' 



(l+t 2 )T 



A plot of the mean current distribution is presented in Figure Q . There are two 
important potential limitations on the applicability of this result: (i) we have 
ignored Coulomb interactions between electrons, which, although they are often 
assumed not to change the single-electron picture of the plateaux transition, may 
well affect the current distribution; (ii) we have neglected quantum tunnelling 
between neighbouring regions of zero potential, which are believed to be relevant 
and to change the universality class away from classical percolation. However, 
for the spin quantum Hall transition|52| this is not the case. 
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Current distribution 
10b 




Figure 7: Average current distribution for the continuum percolation model 

6 Reverse Engineering SLE(k, 2) 

Now let us consider the more general case of two curves starting at the points 
x\ 7^ x 2 - -Pioft, -Pmiddic and Pright are given by the ratio of correlation functions 

(0 2 (a;i)02(a;2)^2(oo)02(oo)) F 1 (x 1 ,x 2 ) 
Using the following definition of the differential operator D\ : 



D 1 



' k d h 2 ,i 1 d 



-Adx\ {x 2 — x\) 2 x 2 ~x\dx 2 
Fq, satisfies the BPZ differential equation 

D x F< s ,(x l ,x 2 ,$) =< (&*&(i)fa(x\)fa(x2)<&(ao)) , 
with S Xl = — — — . 

Writing F$ = P(xi,x 2 ,$)F 1 (xi 7 x 2 ) and using DiF 1 (xi,x 2 ) = 0, 

„ n , -r \ K dFx(xi, x 2 ) dP(x\, x 2 , $) . „, ^> 

D 1 P(x 1 ,x 2 ,$) + —- T a a = 8 Xl P(x u x 2 ,$) . 

2F 1 (xi,X2) oxi ox\ 

^ (35) 

F\(x\,x 2 ) is a three point function. This can be seen by taking the fusion of 
the operators at infinity, ie §\ (oo) = ^3(00). Conformal field theory may be 
employed to fix the form of the three point function up to a constant 

F\(x\,xi) = 



\Xi - X2\*"- 2 ' 1 ~" 3 > 1 \x x ~X 2 \- 2 l K 

From this expression for -Fi(aii, x 2 ), it can be seen that 

dF 1 (x ll x 2 ) _ 2 Fi(xi,x 2 ) 
dx\ k \xi — x 2 \ 



20 



Substituting this into Eq. I|35|) . we obtain 

~ 4% + —^—(4- ~ J-)]g(xi,x 2 ,<S>) = <f xl G(x 1>a : a> $). (36) 
2 otj xi — x 2 0x2 oxi J 

This has the form of the adjoint Fokker-Planck equation corresponding to the 
stochastic process 

1- 2 
dx\ = v ndB t H d< 

xi - a; 2 

2 

cfe 2 = di . (37) 

x 2 - x\ 

This is called SLE(k, 2) [21]. Since the differential equation resulting from this 
choice of stochastic variables is the same as the CFT solution to the 2 curve 
problem, we conjecture that this stochastic process gives the driving term in 
the Loewner equation for one curve, given the existence of the other. Recently 
Dubedat|5H] has argued that such a description follows from the requirement 
that the generators of the Loewner processes for the two curves should commute. 

It is interesting to take k = 0, which results in a deterministic Loewner equation 
with analytic solution. Solving Eq. 1|37[) yields the following forcing function 



x 1 = S j2t+ 5 -^, (38) 

where t parametrises distance along the curve and 5 is the distance between the 
starting points of the curves on the real axis. Kadanoff et al. derived solutions 
to the Loewner equation for various driving terms |27| . including x\ oc yft. The 
solution below is similar to the case of a square root forcing term. The mapping 
from the upper half plane with boundary curves grown up to time t back to the 
upper half plane g t is the solution to the Loewner equation 

^ 1—! < 39 > 

subject to the boundary condition gt=o(z) = z. A change of variables to 
G t = -^=, r = in(£ + t), 



leads to the equation 



dG (G-y+)(G-y.) 



dr " 2(v / 2-G) 
with y + = 2\pl and j/_ = —y/2. In terms of 



(40) 



= 4V2HG + V2) + 2V2HG - 2V2) 
3V2 



the equation may be written 
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with solution 

a S 2 
-H( yt - ) = ln( h t) + const . 

t+f 8 

The constant may be set by the requirement that gt=o{z) = z, then 

-Hi-*) = Hl + %)-H(^). (42) 



The boundary curves z c (t) are the line of singularities which are found by setting 
g t (z) = V2^t + S 2 /8, namely 

-J^-^l + (43, 
After substitution for H(G) this may be simplified to 

4(^) 3 -3^ + /(i)=0, 
with f(t) = 2(1 + 8t/S) 3 / 2 — 1. The solutions are hyperbole of the form 

4a 2 _ ! & 2 = S 2 , (44) 
3 

where the location of the tip is given by z c = a + ib. The limit 5 — > is the 
case which we have quoted in Sec. 14. II In this limit, the curve is a straight line, 
proceeding at an angle of 7r/3 from the positive real axis. By symmetry, the 
other curve is also a straight line, making an angle of 7r/3 from the negative real 
axis. Figure displays the solution. 



7 Summary 

This paper has described how the conformal field theoretic treatment of SLE 
may be generalised to two curves. In the limit that both curves originate from 
the same point, the equations for the probability that a point lies to the left, 
right or between the two curves simplify to a third order ordinary differential 
equation. This is the limit which has been investigated in this paper. Results 
have been obtained for the range < k < 8 in terms of integrals of hyperge- 
ometric functions. The special cases of k = 0,2,4,8/3,8 allow exact analytic 
solutions in terms of elementary functions. 

The application of the result for k = 6 to the quantum Hall problem has been 
explained, along with its limitations. 

It would be interesting to investigate the generalisation of the work in this paper 
to the case of n curves starting from the origin. Although no more difficult in 
principle, the mathematics would be complicated; the solutions are those of 
(n + l)th order ordinary differential equations. 
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A Level 3 null states 



The Virasoro generators have the following commutation relations: 

[L m , L n ] = (m - n)L m+n + -^m(m 2 - l)5 m+n . 

The aim of this section is to find the conditions for an operator to have a null 
state at level 3, which is to say that 

L„(X_ 3 + M^-2i-i + AL^M = . (45) 

Choosing n = 1, 2 in the equation above, leads to equations which determine /x, 
A and h 3 , the scaling dimension of 0, as functions of of the central charge, c. 
First, acting with L\, dropping the |</)) for clarity: 

= LiL_3 + /xiiL_ 2 i-i + \LiL i _ l 
= (4 + 2nh 3 )L_ 2 + (3/J + 6A + &\h 3 )L\ . 

The coefficients of L_2 and L 2 _ 1 must both vanish, since otherwise this would 
imply a null state at level 2. Hence 

4 + 2/i/i 3 = 
3{i + 6A(1 + h 3 ) = . 

These simultaneous equations have solutions: 

2 x 1 



ha' h 3 (l + h 3 )- 

Consider Eq. I|45|l with n — 2 to obtain the dependence of h 3 on the central 
charge, c: 

= L 2 L_ 3 + /j.L 2 L_ 2 L-i + XL 2 L 3 _ 1 
= (5 + y + 4^3 + 4// + I8A/13 + 6A)L_i . 



Hence, 



5 + M(| + 4/i3 + 4) + A(18/i 3 + 6) 



Substituting for /i and A from above, 

5-f (|+4fe 3 + 4)+ 1 (18fr 3 + 6) = 0. 
ft 3 2 /i 3 (l + /i 3 ) 

which is the following quadratic equation in h 3 : 

3hl + h 3 (c -7) + c + 2 = 0. 

This has solution 



7-c± ^/(c- 7) 2 - 12(c + 2) 



(3 



(46) 



Hence, ^, A and h 3 are all restricted to given functions of c, the central charge, 
or equivalently in terms of k, the SLE K variable: 
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B From the correlation function to the differen- 
tial equation 



Consider the correlation function P defined by: 

p = lim ($(Ofo(si)^(y)) _ 
y-°° (<p3(xi)<Pl(y)) 

The contour integral associated with the raising operator L_3 acting on the 
state at x\ can be deformed continuously until it surrounds £. Hence 

= (eL_ 3 ^3(aJi)$(f)^(cx>)> + + 77 ! ^ )^M) • 

(? - zi) 

Using the level three null state condition (see Appendix El, 

(L_ 3 + pL_ 2 L-i + \Li t )fa = , 
this equation may be written as 

= -(enL-tL-Kh*) - (eA£l 1 ^)+e2Re[ ,. 1 , 2 ^1p, (47) 

where the operators at infinity and the operators' dependence on position have 
been dropped for clarity. The first term, involving L-iL-x, can be re- written 

as 

J z — x\ 2iti 

where a cancelling minus sign has appeared from reversing the direction of 
the contour from clockwise to counter-clockwise. The integral is equivalent to 
— > <£(£'), using that the scaling dimension of <I> is zero and defining £' as 

= 



Writing £ in terms of real and imaginary parts as £ = u + i v, 

(u - xi) - v i 
o£ = e- — - . 

[u — X\ Y + v z 

Then, the contribution to Eq. (|47(l is 

d ( u d v d \ 

The contribution from the second term in Eq. I|47|) may be written as 
-eA(L 3 _ 1 3 ^) = P(u - x x ) = ^g^P(^ - *i) > 

where we have used that the function P{x\ 1 u, v) must be a function of P(u — 
x x ,v), ie. 

dP{u~xi) dP(u — xi) 
dxi du 
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Lastly, we consider the contribution from the third term in Eq. I|47f) . Again 
writing £ = u + iv, 



5£ = 



(u — xi + iv) 2 



2uvi 



= €- 



(u 2 — v 2 ) 2 + 4u 2 v 2 (u 2 — v 2 ) 2 + 4u 2 v 2 ' 
where x\ has been set as the origin. The contribution to Eq. 1)47(1 is 



d 



2uv 



8 



(u 2 — v 2 ) 2 + Au 2 v 2 du (u 2 — v 2 ) 2 + Au 2 v 2 dv 



d 



2uv 







(u 2 + v 2 ) 2 du (u 2 + v 2 ) 2 dv 
Putting all three terms together and using P = f(u — x\), 



+ v 2 \ du dv) du du 3 



Next make the following substitutions: 



d \ d d 3 {u 2 -v 2 )-f -2uu#- 



p 



o. 



9P{t) 



9P(t) 



1 dP(t) 
v dt 

t dP(t) 
v dt 



Then the differential equation can be rewritten as 



VfJ, 



(i 



-t d 2 



ufi/v 2 <P 

u 2 + v 2 dt 2 ' u 2 + v 2 V v 2 dt 2 
2uv 



I d 
' ~tf~dt 
u d 



X d 3 
+ ^dt 3 

P(t) = 0. 



(u 2 + v 2 ) 2 v dt (u 2 + v 2 ) 2 v 2 dt 

Multiplying both sides by v 3 the ordinary differential equation becomes 

at d 2 a d , d 3 t 2 -1 d 2t 2 d 

+ A- 



- 2- 



t 2 + ldt 2 t 2 + ldt dt 3 {t 2 + l) 2 dt (t 2 + l) 2 dt 
Collecting terms, this is: 

d 3 P 2[xt d 2 P (3 - pt)t 2 - (1 + fi) dP 



P(t) = 0. 



X 'dt 3 t 2 + 1 dt 2 
which is Eq. (|f 3(1 in the text. 



(t 2 + iy 



dt 



0, 
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